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We report a de Haas-van Alphen (dHvA) oscillation study of the 111 iron pnictide superconductors LiFeAs 
with Tc ~ 18 K and LiFeP with Tc ^ 5 K. We find that for both compounds the Fermi surface topology is in 
good agreement with density functional band- structure calculations and shows quasi-nested electron and hole 
bands. The effective masses generally show significant enhancement, up to ^ 3 for LiFeP and ^ 5 for LiFeAs. 
However, one hole Fermi surface in LiFeP shows a very small enhancement, as compared with its other sheets. 
This difference probably results from fc-dependent coupling to spin fluctuations and may be the origin of the 
different nodal and nodeless superconducting gap structures in LiFeP and LiFeAs respectively. 



Identification of the particular structural and electronic 
characteristics that drive superconductivity in the iron-based 
materials continues to be a central experimental and theoreti- 
cal question in the field. A successful theory needs to explain 
trends, such as the variation of Tc and also the structure of the 
superconducting energy gap. In most of the iron arsenides the 
parent materials have a non-superconducting, antiferromag- 
netically ordered ground state. Disruption of this magnetic 
order leads to superconductivity and then eventually a non- 
superconducting paramagnetic ground state. A good exam- 
ple is the BaFe2(Asi_a;Pa;)2 series which has a maximum Tc 
=30 K when x o:^ 0.33 dill. Here BaFe2As2 has a magnetic 
ground state whereas BaFe2P2 is a paramagnet and neither 
superconduct. 

The Ill-family of iron-pnictides LiFeAs i-a^P^^, is unique 
because both LiFeAs and its counterpart LiFeP superconduct 
and are non-magnetic with Tc ~ 18 K f3l,'4l and ~ 5 K |5|, re- 
spectively. Also, penetration depth measurements have shown 
that LiFeAs is fully gapped (61 [3, whereas LiFeP has gap 
nodes Q. Establishing whether this switch of pairing struc- 
ture is linked to changes in the topology and orbital character 
of the Fermi surface (FS) provides an stringent test of candi- 
date theories for the superconducting pairing in these materi- 
als. 

Magneto-quantum oscillation effects such as the de Haas- 
van Alphen (dHvA) effect are a powerful probe of the three- 
dimensional bulk Fermi surface and have been successfully 
used to study a variety of iron-based superconductors IS] O . 
In this Letter, we present a study of the dHvA oscillations in 
both LiFeP and LiFeAs which establishes that the bulk Fermi 
surface topology of these compounds is in good agreement 
with DFT calculations. Furthermore, by comparing the val- 
ues of the extracted effective masses of the quasiparticles to 
the calculated band masses, we find significant orbit depen- 
dence to the mass enhancement factors which is likely linked 
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FIG. 1. (color online) Torque versus field for LiFeP and LiFeAs. 
The top panels show the raw pulsed field torque data in units of the 
change in the cantilever resistance at T = 1.5 K. The arrow indi- 
cates the position of the irreversible field. The middle panels show 
the oscillatory part of the torque after subtraction of a smooth back- 
ground. The bottom panels show FFTs of the torque. For the peak 
labels see the main text. For LiFeP we show FFT spectra computed 
over different field windows (a) (25-58 T) which shows the splitting 
of the a peaks, (b) (40-58 T) which decreases the influence of noise 
on the higher frequency peaks, and (c) (33-45 T) for the do field data 
at T = 0.45 K and = 51°, showing the strong ( peak. 



to the contrasting superconducting gap structures and Tc in 
these compounds. 

Single crystals of LiFeP and LiFeAs were grown by a flux 
method 1 10|. Small single crystals, typically 50 x 50 x 10 /im^ 
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FIG. 2. (color online) (a) Calculated Fermi surfaces of LiFeP. (b) and (c) Show the evolution of de Haas-van Alphen frequencies with magnetic 
field angle. Experimental data are shown in the right panels as symbols (triangle = pulsed field, square/circle=dc field; circles=probably 2nd 
harmonics). The solid lines show the result of the DFT calculations; the bands are shifted in the right hand panels to best fit the experimental 
results. The numbers refer to the bands in (a). In all panels the frequencies have been multiplied by cos for clarity.(d) Slices through the 
determined Fermi surfaces at particular kz values (with shifted bands). The dashed/solid lines are the hole/electron sheets respectively, and the 
latter have been shifted along the [110] direction such that their center coincides with the holes. 



for LiFeP and 100 x 200 x 50 /im^ for LiFeAs, were selected 
for the torque measurements. To avoid reaction with air the 
samples were encapsulated in degassed Apiezon-N grease. 
Sharp superconducting transitions were measured using ra- 
dio frequency susceptibility with Tc onset (midpoint) values 
of 4.9 K(4.7 K) and 18.4 K(17.3 K) for LiFeP and LiFeAs, 
respectively. The samples were mounted onto miniature Seiko 
piezo-resistive cantilevers which were installed on a rotating 
platform, immersed in liquid ^He, in the bore of a pulsed mag- 
net up to 58 T in Toulouse. Measurements on the same crys- 
tals were also conducted in an 1 8 T superconducting magnet 
in Bristol and a 33T Bitter magnet at HMFL in Nijmegen 
and 45 T hybrid magnet at NHMFL, Tallahassee, all equipped 
with ^He refrigerators. 

Torque versus magnetic field data are shown in Fig. [T] For 
both materials dHvA oscillations are seen at high fields and 
low temperatures, well above the upper critical field, esti- 
mated to be < 1 T for LiFeP El and -16 T for LiFeAs ifHll 
when B\\c (see also Fig. [T]). After the fast Fourier transform 
(FFT) as a function of inverse field, several strong peaks are 
visible (Fig. [T] bottom panels), which correspond to the ex- 
tremal cross-sectional areas A^ of the FS: F = hAk/27Te. 
For LiFeP, the spectrum is dominated by two low dHvA fre- 
quencies around 300 T and 400 T, labelled ai and 0^2. The 
amplitude and frequency of the peak at '^ 750 T is consistent 
with this being the second harmonic of the a peaks. The other 
five peaks (/3, 7, (5, £, are clearly derived from unique Fermi 
surface orbits. For LiFeAs, three frequencies are visible at 
1.5 kT, 2.4 kT and 2.8 kT, labelled as /3, 5, s, respectively. 

To properly identify these FS orbits, we performed field 
sweeps with different field orientations starting from 6> = 0° 
(B\\c) and rotating towards the a6-plane. For a perfectly two 



dimensional (2D) FS, F ex 1/cos^, so by multiplying F by 
cos 6 the degree of two dimensionality of a FS can be easily 
seen. For quasi-2D surfaces, F cos will decrease with in- 
creasing for a local maximum of Fermi surface orbit area as 
a function of kz whereas the opposite will be true for a local 
minimum. The data in Fig. |2] suggest that for LiFeP e and 7 
are a maxima, and (3 and S are minima. The two lowest fre- 
quency a orbits have opposite curvature indicating that they 
are the maximum and minimum of the same FS sheet. At an- 
gles close to 6> = 50° strong peaks are seen (labelled ( in Figs. 
[T]and[2]) which are likely from the outer hole sheet (band 3). 
The amplitude becomes large at this angle because of the Ya- 
maji effect, expected when the two extremal orbits of a quasi 
two dimensional Fermi surface cross. Close to = 30°, an 
additional branch r] is visible. For LiFeAs, the e orbit is a 
maximum, while S and (3 orbits are likely to be minima orbits. 

To identify the origin of the orbits and solve the structure 
of the Fermi surface we have performed DFT calculations us- 
ing the linear augmented plane wave method, implemented 
in the wien2k package 1 13 1. We used the experimental crys- 
tal structure 1 14 1 and included spin-orbit coupling (SOC). The 
calculated Fermi surfaces (see Fig.|2ja,d)) are quite similar for 
both materials, there are three hole bands at F and two elec- 
tron bands at M as found previously {TSJ . The two outermost 
hole sheets are quite 2D, whereas the innermost xz/yz hole 
pocket is strongly hybridized with ^^2 near Z and is closed 
there, while remaining 2D away from this point. By con- 
trast, the electron orbits are very strongly warped. This ge- 
ometry is reflected in the calculated angular dependence of 
the dHvA orbits (Fig.|2];b,c)). For the 2D hole sheets F cos l9 
varies little with angle and the maximal and minimal area 
are close. For the electron sheets there is a large deviation 



from this behavior. For LiFeP, SOC splits the two outer- 
most hole bands, which are accidentally nearly degenerate in 
non-relativistic calculations, and causes their character to be 
mixed dj-z/yz/dxy In LiFeAs these bands are well separated 
irrespective of SOC and have a predominantly d^z/yz (mid- 
dle) and dxy (outermost) character. The SOC also splits the 
electron bands along the zone edge (X-M) inducing a gap 
of ~35meV (see [2jd)), hence as in LaFePO |16| we esti- 
mate that magnetic breakdown orbits, along the elliptic elec- 
tron surfaces in the unfolded Brillouin zone, to be strongly 
damped. 

By comparing the calculations to the data (Fig. |2jb,c)), in 
particular the curvature of F cos 0, the correspondence be- 
tween the observed dHvA frequencies and the predicted Fermi 
surface orbits is immediately apparent for most orbits. The 
observed (3 frequencies are likely a mixture of signals from 
orbits 2a (hole) and 5a (electron) close to 6> = 0° but are sep- 
arately resolved at angles close to 30° (the r] branch probably 
corresponds to band 5a). For LiFeP, relatively small shifts 
(somewhat smaller than for other Fe pnictides |9|) of the band 
energies: +20 meV and +45 meV for band 4 and 5 (electron) 
and —65, —80, 18 meV for bands 1, 2 and 3 (hole) bring the 
observations and calculations into almost perfect agreement 
as shown in Fig.[2jb). As in other Fe pnictides ||9j|T2l, these 
shifts shrink both the electron and hole FSs and likely orig- 
inate from many body corrections to the DFT bandstructure 
ifTSJI . Although the maximal orbit of band 4 which is close to 
6 kT was not observed, probably because the scattering rate in 
our sample was too high, we can estimate the accuracy of our 
band energy determinations by calculating the difference in 
total volume of the electron and hole Fermi surfaces. We find 
a small imbalance of just +0.02 electrons per unit cell which 
shows the consistency of the procedure. 

For LiFeAs, the curvature and absolute values of Fcos 
suggest that the e orbit originates from the maximum of the in- 
ner electron Fermi surface (band 5) and the extended angular 
dependence of the S orbit suggest that this originates from the 
minima of the electron surface (band 4a), rather than the max- 
imum of the middle hole surface (band 2b) which is of similar 
size in the calculation. The limited angular extent of the data 
for the /3 orbit means it is not possible to say if it originates 
from band 5 a (electron) or band 2a (hole) although 5a seems 
more likely. To accurately match the e and S orbits with the 
calculations only very small shifts of the band energies are re- 
quired (—5 meV and +18 meV for bands 4 and 5 respectively) 
(Fig. |2jb)). We did not observe the smallest hole FS (band 
l)in LiFeAs, even though the same band gave the largest sig- 
nals for LiFeP. This suggests that band 1 does not cross the 
Fermi level in LiFeAs, which requires that it shifts down by 
^ —40 meV, possibly because of enhanced SOC. The small 
shift of the electron bands is almost perfectly compensated by 
the removal of band 1 , so the remaining hole bands are not 
shifted in Fig. |2jb). Although we do not see definitive evi- 
dence for the hole orbits, probably because of a significantly 
higher impurity scattering rate in LiFeAs compared to LiFeP, 
the small size of the energy shifts needed to match the elec- 
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FIG. 3. (color online) Quasiparticle effective masses determination. 
Amplitude of the FFT peaks (the field ranges as indicated) versus 
temperature. The lines are fits to the Lifshitz-Kosevich formula I.20J . 
The effective mass values are shown in Table U 



tron bands combined with the similar small shifts required in 
LiFeP to match both electron and hole bands strongly suggests 
that the DFT calculations correctly predict the Fermi surface 
topology of these 111 compounds. This is contrast to the pho- 
toemission results of Borisenko et al 1 19] for LiFeAs, where 
a significant discrepancy between the size of the hole sheets 
and the DFT calculations was found. 

The strength of the electron-electron interactions can be es- 
timated from measurements of the quasiparticle effective mass 
m* on each orbit through the temperature dependence of the 
amplitude of the dHvA signals, by fitting the latter to the 
Lifshitz-Kosevich formula 1201 (Fig.|3]). These measurements 
were conducted in dc field on the same samples to avoid any 
possibility of sample heating at low temperature. The derived 
values along with the DFT calculations are shown in Table [l| 

For LiFeP, the enhancements factors A = m*/m5 — 1 vary 
strongly between orbits. For the electron sheet A is in the 
range 1.4-2.3, which is comparable to values found for the 
electron sheets of LaFePO (Tc=6K) II2TII . The smallest and 
largest hole sheets (bands 1 and 3) are also strongly enhanced, 
however for the middle hole sheet (orbits 7 and /3, band 2) A 
is ~ 3 times smaller than for the other sheets, despite having 
similar orbital character. As an enhancement A :^ 0.2 |22J is 
expected from electron-phonon coupling, this means that the 
residual electron-electron component for this particular orbit 
is very weak. This is an interesting observation, relevant to 
the ongoing discussion ll23ll as to whether the mass enhance- 
ment comes entirely from local correlations or partially from 
long range spin fluctuations. If the mass renormalization in 
this compound is due to the same spin fluctuations that are be- 
lieved to cause superconductivity, we can conclude that band 
2 is very weakly coupled with these fluctuations, so that the 
pairing amplitude on this band will be small and hence it is a 
possible candidate for the location of the gap nodes. Calcu- 
lations suggest 1 23 1 that node formation is controlled by the 
xy pocket, so that if this pocket exists, the order parameter is 
nodeless, otherwise nodes form on an electron (band 4, in our 
notation) pocket. LiFeP seems to deviate from this rule, as it 
has a well developed xy pocket (band 3) and has gap nodes. 
LiFeP therefore appears to be a challenging and an extremely 
interesting material for further theoretical modelling. 



TABLE L Measured and calculated dHvA frequencies. The mea- 
sured frequencies are extrapolated to ^ = 0. The effective (m*) and 
calculated (mb) band masses are quoted in units of the free electron 
mass (— sign indicates hole orbit). For the LiFeP the values marked 
with t were determined at an angle of 51 ° (A was calculated with 
nib also calculated at this angle). At this angle orbits 3a and 3b cross 
(and have maximum amplitude), and orbit 2a is clearly differentiated 
from 5a. Values marked with [j potentially have overlapping orbits 
and so the mass enhancements and orbit assignments are less certain. 
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For LiFeAs, the measured effective masses are uniformly 
larger than in LiFeP. For the electron sheet (band 5) A is more 
than 3 times larger than in LiFeP. This observation suggest 
that mass renormalization in iron pnictides is caused by the 
same interaction that drives superconductivity. This agrees 
with previous findings in the isoelectronic superconducting 
series, BaFe2(Asi_a^Pa^)2, in which the effective mass of the 
electron bands are closely related to the increase in Tc ifTTl . 
Interestingly, the large mass enhancement in LiFeAs is not 
accompanied by a corresponding large shrinking of the Fermi 
surface volume ifTTl . 

In summary, dHvA oscillations have been observed in two 
members of the 111 family of superconductors, LiFeP and 
LiFeAs. In both cases we find that measured data are con- 
sistent with the topology of the DFT calculated Fermi surface 
with small band energy shifts. The many-body mass enhance- 



ments are larger in LiFeAs than in LiFeP. In LiFeP, the mid- 
dle hole band has significantly lower mass enhancement than 
the other sheets, which implies that the electron-hole scatter 
rate is suppress for this sheet. This may be the origin of the 
lower Tc and nodal gap in LiFeP, and suggests that the mass 
enhancement is to a large extent due to a /^-dependent spin- 
fluctuation induced interaction, which are also responsible for 
the pairing. It will be very interesting to see whether these fea- 
tures and the nodal gap structure in LiFeP can be explained by 
detailed microscopic calculations. 
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